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Chapter 1 



The Mean Field Theory of Nuclei 
and The Bag Model of Baryon 



1.1 Mean Field Theory of Nuclei (MFT) 

Two model field theories of the nuclear system were studied in detail by Serot and Walecka 
[1]. The first is based on baryons and neutral scalar and vector mesons (model QHD-I). 
The quanta of the fields are the nucleons (n,p) and the sigma (a) and omega (cu) mesons. 
The neutral scalar meson couples to the scalar density of baryons through g s ipip(() and the 
neutral vector meson couples to the conserved baryon current through gvip^^ipV^. 

The Lagrangian density for this model is 

C = ^[-fiidn-gvVj-iM-gsfiW 

+\ (d^m - m 2 s 4> 2 ) - + \m 2 v V»V», (1.1) 

where 

= d/jV,, - d„V,. 
Lagrange's equations yield the field equations: 

[^{id,-g v V,)-{M-gM^ = 0,_ (1.2) 
{d^-m 2 s )<t> = g s H>, (1.3) 
d^-mlV = g v W1>, (1-4) 

Eq. (1.2) is the Dirac equation with the scalar and vector fields. Eq. (1.3) is simply 
the Klein-Gordon equation with a scalar source. And Eq. (1.4) looks like massive QED 
with the conserved baryon current, 

A u =^Y*P- (1-5) 
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The energy-momentum tensor is: 



2 



-<9 A 0«9 A + m 2 s( p 2 + l -F Xu F x ° - m 2 v V x V x 



+iip^d u ip + d^(j)d u (j) + d u V x F Xtl . (1.6) 

Lagrange's equations ensure that this tensor is conserved and satisfies d^T^ = 0. It 
follows that the energy-momentum P u defined by 



P v = I d'xT^ (1.7) 



is a constant of motion. 



We observe that at high baryon density, the scalar and vector field operators can be 
replaced by their expectation values, which then serve as classical, condensed fields in 
which the baryons move, 

- (0)=0 O , (1-8) 
V, - (V^=S^V . (1.9) 

For a static, uniform system, the quantities 0o an d Vq are constants independent of 
x^. Rotational invariance implies that the expectation value (V) vanishes. 

The condensed, constant, classical fields 0o and Vq are directly related to the baryon 
sources. The source for Vq is simply the baryon density p A = A/V. Since the baryon 
current is conserved, the baryon number 

A= d 3 x^ifj (1.10) 
Jv 

is a constant of motion. For a uniform system of A baryons in a volume V, the baryon 
density is also a constant of motion. In contrast, the source for 0o involves the expectation 
value of the Lorentz scalar density {ipijj) = 0o- This quantity is dynamical that must be 
calculated self-consistently using the thermodynamic argument that an isolated system 
at fixed A and V (and zero temperature) will minimize its energy, 

dE(A, V; O ) 

1 ' ; =o. in) 

d<p 

In the MFT, the Lagrangian density is 

%FT = ^ [i-l^d^ — 9vl°V — (M — g s (/>o)] ip 

-\<<Pl + \rnlVl (1.12) 

Hence, the Dirac equation is linear, 

[i^O, - g vl °V — (M — g.^)} ^ = 0, (1.13) 
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and may be solved directly. We seek normal- mode solutions of the form ip(x^) = ip(x)e lEt . 
This leads to 

Hip(x) = Eip(x), (1.14) 
H = [-iaV + g v V + f3(M-g s (f) ). 

The effective mass M* is defined by 

M* = M-g s <j) . (1.15) 

The condensed scalar field 0o thus serves to shift the mass of the baryons. Evidently, the 
condensed vector field V shifts the frequency (or energy) of the solutions, 

E = g v V + E*. (1.16) 

In the MFT, the energy-momentum tensor is 

It then follows that 

P"= f d?xT 0v . (1.18) 
Jv 

The second nuclear model is more realistic. To discuss nuclei with Z ^ N, it is 
necessary to extend model QHD-I to include p mesons, which couple to the isovector 
current, and the coulomb interaction. 

Since we are also interested in comparing quantitative predictions with experiment, we 
will extend the model to include the p meson and photon fields (QHD-II). The Lagrangian 
density is 



C = i> 



7" (id? - 9vV^ - ^gpr.b^ - X -e{\ + r 3 )A^j — (M — g s <p) 



1> 



-\g, v .& v + \m%& - \H^\ (1.19) 



where 



F 


= dpVv 


- d v V„, 


—* 


= dj) u - 


- d^bp, 




= d,A v 


- d v A^ 



The field equations (1.2 to (1.4) must also be extended to include contributions from 
the rho and photon fields, 

(d^-m% = grfn^, (1.20) 
UA, = e^l±^ 7 ^.21) 



7 M (id„ - g v V„ - l -g p r\ + T*)A^j — (M — g s cf)) 



i\) = (1.22) 
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If N Z, the neutral field p° corresponding to can develop a classical, constant 
ground-state expectation value in nuclear matter according to 

(b^) = 5^5 j %. (1.23) 
The baryon field, however, obeys a Dirac equation analogous to (1.13) namely, 

il^ - g vl V - \g p T zl % - l -e(l + r 3 ) 7 % — (M — (? S O ) V> = (1.24) 



Although the baryon field is still an operator, the meson fields are classical; hence, Eq. 
(1.24) is linear, and we may seek normal-mode solutions of the form %p{x^) = ip(x)e~ lEt . 
This leads to 



hip{x) = Eip(x) 
h = 



1.25) 



ia.V + g v V + \g p T Z b + \e{\ + r 3 )A + (3 (M — g s <j) ) 



which defines the single-particle Dirac Hamiltonian h. 

The single-particle wave functions in a central, parity-conserving field may be written 



as 



1p a (x) = i)nkmr{x) 



km 



iy> rullti 



Cr 



;i.26) 



The equations for the baryon wave functions follow immediately upon substituting 
1.26) into (1.25): 



dG a k 

—j 1 — 

ar r 
dF a k 

~j <J"a + 

ar r 



E a - g v V - ^g P r a b - ^e(l + r a )A + M - g s (p 
E a - g v V - ^g P T a b - ^e(l + r a )A - M + g s (p 



0, (1.27) 



G a = 0. (1.28) 



For a given set of meson fields, the Dirac equations (1.27) and (1.28) may be solved 
by integrating outward from the origin and inward from large r, matching the solutions 
at some intermediate radius to determine the eigenvalues E a . Analytic solutions in the 
regions of large and small r allow the proper boundary conditions to be imposed. 



1.2 Relativist ic Nuclear Physics 

Relativistic nuclear physics was first pioneered and developed by Shakin and Celenza [2]. 
The systematic application of this theory is able to resolve some long-standing puzzles in 
the theory of nuclear structure such as the binding energy and the saturation density of 
nuclear matter, the effective force in nuclei, and the nucleon self-energy for bound and 
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continuums. Its success is based upon the use of the Dirac equation for the description 
of motion of a nucleon. The potentials appearing in the Dirac equation are assumed to 
contain large (Lorentz) scalar and vector fields. 

The scalar fields enter the Dirac equation in the same way as the nucleon mass. Since 
these fields are quite large (—400 MeV), they have the effect of including a major reduction 
of the nucleon mass when the nucleon is in the nuclear matter. It is the description of this 
change of mass of the nucleon that is an essential element in the success of the relativistic 
approach. It is further necessary to understand that the vector field seen by a nucleon 
is large and repulsive, so that the energy of the nucleon in the nuclear matter does not 
differ very much from the energy of a nucleon moving in the weak fields which appear 
in the standard Schrodinger description. More precisely, the description relation relating 
the energy and momentum of a nuclear quasiparticle is similar to that of the Schrodinger 
theory. Therefore, the system may be said to exhibit "hidden" relativity. 

Indeed, for decades the Schrodinger approach to nuclear structure physics provided a 
reasonably satisfactory model of nuclear dynamics. It is only in the last decade that the 
true relativistic features of the system have become apparent. 

We write the Dirac equation for a nucleon in the nuclear matter as 

[a.p + 7°m + V(p)](j)(p, s) = e0(p, s), 

where V(p) is the potential. It will be useful to introduce the self-energy = 7 t V(p) 
and rewrite this equation as 

[f.p + m + S(p)]0(p, s) = 7°e0(p, s). 



Now let us assume that the self-energy is of the form 
so that we have 

[l-P + (m + AMp, s) = 7 °(6 - B)cj>(p, s). 
A positive-energy spinor solution of this equation is 

\E{p) J 

where 

E(p) +fh ( Xj 

2 ™ V wl^ Xs 



u(p, s, m) 



Here u(p, s, rh) is the positive-energy solution of the Dirac equation without interac- 
tion, except for the fact that the nucleon mass m has been replaced by m = m + A and 
E(p) = ^/p 1 + rh 2 . The normalization chosen here is 



so that 

0V = i- 

We further note that the energy eigenvalue is 

e = B + + m 2 

-2 

= m + 5 + A + ^ + ---. 

2m 

Now, as we have mentioned, A is large and negative (—400 MeV) and B is large 
and positive (300 MeV). Therefore A and B largely cancel and the dispersion relation 
is essentially the same as that which one could find in a nonrelativistic model. The 
development of this simple relativistic nuclear model leads to two categories. The first 
we will call Dirac Phenomenology. This category is distinguished by having several free 
parameters which are adjusted to fit nuclear date. The second category will be called 
Relativistic Brueckner-Hartree-Fock (RBHF) theory and is characterized as having no 
free parameters other than those introduced in fitting free - space nucleon - nucleon 
scattering data. Interest in the development of the RBHF approximation grew out of 
the successful application of Dirac phenomenology to the description of nucleon - nucleus 
scattering data. 

We have so far presented only the main idea and the materials necessary to our con- 
sideration. For those who are interested to the detailed results of the relativistic nuclear 
physics, please, see the monograph of Celenza and Shakin quoted above and the references 
herein. 



1.3 Bag Model of Baryons 

It is well accepted QCD is the theory of strong interactions. However, in general, QCD is 
never solvable: at low energies and small momenta transfer the running coupling constant 
a s > 1. The bag model, outlined for the first time by the group of M.I.T. theorists [3], 
is a phenomenological approach, in which two basic features of QCD are incorporated: 
asymptotic freedom and confinement. 

The main assumption of the M.I.T. bag model states that, baryon is considered to be 
a bag of spherical shape, in which the constituent quarks move freely and are described 
by the Dirac equation 

rr i - dl l> 

H * = l w 

with the Hamiltonian 

H = a.p + (3M. 

Consider the case k — — 1, which is the 51/2 level. The solution of this equation has 
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the form 

/ (-> + \ AT I V ~E~ 3° ill) \ -,tl 

If we parametrize the energy levels as 

E nk = u nk /R, E nk = VE* - M 2 , 
the density of quarks is readily calculated as 



where 



■ 2 (ur\ E-M 2 /c^rx 



i?/ E + M J1 \R) 

1 r < i? 
r > R. 



9 V , 



Thus, the density certainly does not vanish at r — R. Clearly, although the lower 
component is suppressed for small r, it does make a sizeable contribution near the surface 
of the bag. Of course it is natural to ask whether this is not unusual in comparison with 
nonrelativistic experience, where ip(R) would be zero. However, such a solution would 
not be consistent with the linear Dirac equation. What counts is that there should be no 
current flow through the surface of the confining region. For example, in the MIT bag 
model it is required that 

n^^ip = 

at the surface - where n M is a unit four vector normal to the surface of the confining region. 
In the MIT bag model this condition is imposed through a linear boundary condition 



at the surface. This implies 

and hence 

because 



ip = —iip'-y.n, 



Consider now the normal flow of current through the bag surface: 

= —ipip = vpip = 0. 
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Thus, it is not the density, but ipip which should vanish at the boundary in the rela- 
tivistic theory, 

- , E + M 9 . . E-M . 

H\r=R = — e~ ilM ~ — e~ = °" 

That is, the matching condition is exactly equivalent to the linear boundary condition 
(l.b.c.) for the static spherical MIT bag, 

i^.nip = —i^.fip = if), 

where 

rf = (0,f). 

We consider the energy-momentum tensor for a model, 

Ty = T^Oy, 

and T^ v is the familiar energy-momentum tensor for a free Dirac field 

The condition for overall energy and momentum conservation is that the divergence 
of the energy-momentum tensor should vanish, and this is certainly true for T 7 ^, as is 
easily proven from the free Dirac equation 

d^T^ = 0. 

However, the fact that these quarks move freely only inside the restricted region of 
space V leads to problems. Indeed, 

d^v = n^A s , 

where A s is a surface delta function 

A s = -n.d(9v). 

In the static spherical case we find that A s is simply 5(r — R). Putting all these together 
we obtain 

d,?? = i^.nd^As, 

and using the l.b.c. 



d»T{? = -- d u (ipijj)\ s A s = -Pn v A s , 



1 

2 

where P is the pressure exerted on the bag wall by the contained Dirac gas 

P = - l -n.d»{W)\ s . 

Clearly, this model violates energy-momentum conservation. Furthermore, this violation 
is an essential result of the confinement process. 
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The resolution of this problem, we add an energy density term B6y to the Lagrangian 
density. Then (since T^ v involves Cg^) the new energy-momentum tensor ^j^IT ^ as ^ e 
form 

T MIT = ^ + B 9n°v- 
Therefore, the divergence of the energy-momentum tensor is 

d»T^ 1T = (-P + B)n»A s , 

which will vanish if 

B = P = --±n.P(frl;)\ a - 

Therefrom a relativistic bag model of nuclei will be proposed, the A-dependence of 
nuclear radius will be calculated and the Weizssacker formula for nuclear binding energy 
will be derived exactly if the corresponding parameters of the model are adequately fitted. 
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Chapter 2 

Bag Model of Nuclei (Z = N) 



2.1 Main Assumption 

We consider the simplest possible case of A baryons moving inside a spherical volume of 
radius R, outside of which there is a pressure exerted on the nuclear surface. 

Let us therefore begin with the Dirac equation for a baryon of mass M: 

(i 7 "fy - M - E)^) = 0, (2.1) 
where £ is the baryon self-energy having the form 

£ = + 7^%. (2.2) 
Inserting (2.2) into (2.1) we obtain the equation 

[7"(i0„ - Vy - (M + 0)]^) = 0. (2.3) 

Eq. (2.3) is nonlinear quantum field equation and its exact solution is very com- 
plicated. We have therefore made little progress by writing down this equation with a 
suitable method for solving it. 

In the MFT, 

- (0)=0o, (2.4) 
V, - {V„) = S^Vq. (2.5) 

For a static, uniform system the quantities 0o and Vo are constants. Hence, the Dirac 
equation is linear, 

[^-9 M - 7 V o -(M + o )^(x M ) = O, (2.6) 

and may be solved directly. 

Our basic assumption is formulated: the nucleus A is considered to be a MIT bag, 
inside of which the motion of nucleon is described by the Dirac equation (2.6). The 
quantities O and V can be determined only after we have fitted to experimental data. 



12 



Let us next consider the energy-momentum conservation for nuclear bag. For stable 
nuclei, there should be no current flow through the surface of the confining region. In the 
MIT bag model it is required that 

in^ = = ^ifj = (2.7) 

at the surface, where n M is a unit four vector normal to the surface of the confining region. 
Thus, ipip which should vanish at the boundary in a relativistic theory. The matching 
condition of the present model is exactly equivalent to the linear boundary condition for 
the static spherical MIT bag. 

The Lagrangian density for the present model is 

C = #7^ - 7 °Vo - (M + o )]^ y + B9 S , (2.8) 

where B9$ is a energy density term. 

Then the energy-momentum tensor has the form 

T Bag = ^9^6 V + Bg^Os, (2.9) 

where 9 V and 9s define the bag volume and the surface 



0rH J , V ={ 1 r=H (2.10) 



r < R 
r = R 
r > R 



and B is the constant surface tension. 

Therefore, the divergence of the energy-momentum tensor is 

'1 



v Bag 



-n.d(^)\ s + 2B 



= (-P s + 2B)n»A s . (2.11) 
The condition for energy and momentum conservation is 

^ ag = 0, (2.12) 

hence 

B = -±n.d(W)\ s , (2.13) 

where P$ is the pressure exerted on the bag wall by the contained A baryons. 

E. (2.13) involves the square of the baryon fields, and is referred to as the nonlin- 
ear boundary condition of the MIT bag model of nuclei. Because of this condition the 
introduction of a constant surface tension B involves no new parameters. 

Eq. (2.6) and the condition (2.13) constitute the basic ingredients of our model. 
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2.2 Wave Function for Baryon 



The Dirac equation for the present model is linear and may be solved directly. We seek 
normal-mode solutions of the form 

^(x fl )=^(r)e- lEt 6 v . (2.14) 

The Dirac equation then becomes 

Hip(r) = Eip(r), (2.15) 
H = H5.V + K, + /3(M + 0o)]- 

The effective mass M* is defined by 

M* = M + (j) . (2.16) 

The scalar field 0o thus serves to shift the mass of the baryons. Evidently, the vector field 
Vq shifts the frequency (or energy) of the baryon, 

E* = E-V . (2.17) 

Hence, Eq. (2.15) becomes 

(-ia.V + (3M*)ijj(r) = E*^(r). (2.18) 

The single-particle wave functions in a central, parity-conserving field may be written 



as 

1pa(r) = ^nfemr(r) = I k r . (2.19) 



^km 

■ r nkr V) 

' - ^ —km 



; F nkr (r) , 
r 

Their angular momentum and spin parts are simply spin spherical harmonics 

1 1 

$fcm = (^ m «9 m ^-jm)^ mi Xm s , (2.20) 



mi,m s 



k-l l = +U + 1)> ° (2 21) 

\ -(/ + i) = -(j + i)<o ' ^ 21) 

where Yi mi is a spherical harmonic and Xm s is a two-component Pauli spinor. The label 
a, {a} = {a;m} = {n/cr;m}, specifies the full set of quantum numbers describing the 
single-particle solutions. Since the system is assumed spherically symmetric and parity 
conserving, a contains the usual angular-momentum and parity quantum numbers. ( T is 
a two-component isospinor. The principal quantum number is denoted by n. The phase 
choice in (2.19) leads to real bound-state wave functions G and F for real potentials in 
Hamiltonian (2.15). 
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The equations for the baryon wave functions follow immediately upon substituting 
(2.19) into (2.18)0: 

4- + - ] G - (E* + M*)F = 0, (2.22) 
dr r J 

A. - -\ p + (E* - M*)G = 0. (2.23) 
dr r J 

These equations contain all information about the static ground-state nucleus. They 

are coupled linear differential equations that may be solved exactly for a given set of 
potentials. 

Consider the case k — — 1 which is the level Si/ 2 . Eq. (2.22) implies 

F= (E* + M*y 1 ( —-- )G, (2.24) 
\dr r J 

so that defining 

W 2 = E* 2 - M* 2 , (2.25) 
the equation for the upper component of ip a (f) is 

^2 + W2 ) G = ( 2 - 26 ) 
The solution of this equation has the form 

G(r)=C sin Wr, (2.27) 

and hence [from Eq. (2.24)] 

F(r) = C(E* + M*)~\W cosWr - sin Wr/r). (2.28) 

The solutions of the Dirac equation (2.18) come from (2.19) is written as 

Mrl=c( . )(r. (2.29) 

The normalization condition that yields the numbers of baryons contained in the nucleus 

A, 

r d 3 x^ip = A. (2.30) 



Now let us assume that the bag has a spherical shape with radius R. ipif) which should 
vanish at the boundary in a relativistic theory. Eq. (2.29) implies that [see Eq. (2.7)], 



W 2 



Mr=n = fo(WR) - - Mm)2 jf(WR) = (2.31) 



1 We use t?.V(G$fc m /r) = —(l/r)(d/dr + fc/r)G$_fc m , and a similar relation for F. 
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and hence 



J^'^=\/^]^fi( WR )- ( 2 - 32 ) 



This is appropriate boundary condition for confined baryons. Thus, the boundary 
condition of the MIT bag model is used, which provided the eigenfrequency of baryon u a , 
if we parametrize the energy levels (wavenumber) as 

W a = u a /R; {a} = {nkr}, k = -1, (2.33) 

where n is the principal quantum number and u a satisfies the equation [from Eq. (2.32)] 

tancj a = —. =. (2.34) 

1 - M*R - ^u 2 a + M* 2 R 2 



Hence, the eigenvalues E a may be determined by matching the solutions at some 
intermediate radius. Analytic solutions in the restricted region of space V allow the 
proper boundary conditions to be imposed. Taking into consideration (2.33) we get the 
energy spectra for baryon, 



E a = ±^W[TM^+V 



±)l^ + M * 2 + V o- ( 2 - 35 ) 

For convenience, the sign (— ) drops out in what follows. The single-particle wave 
functions now has the form 



jo (f r) *- 



Mr) = C( ff. v *://r^g )( T . (2.36) 

Given the general form of the solutions in (2.36), we may now evaluate the local baryon 
density. Assume that the nuclear ground state consists of filled shells up to some value of 
n and k. This is consistent with spherical symmetry and is appropriate for magic nuclei. 

With these assumptions, the local density of baryons is readily calculated as 



pA = ip "ip9 



•^0 ( 73 / 7-1* i Jl/f* 



R J E* + M* \ R 



(2.37) 



Substituting Eq. (2.37) into Eq. (2.30), we can calculate the normalization constant 
which is defined by Eq. (2.30) for k = —1, 

C 2 = ^ ( E* a + M* \ E* a (E* a -M*)R 

" 4nR*j$(u a ) V E* J 2E?R - 2E* + M* ' { ' ] 



2 We use (2.25) and (2.33). 
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Finally, the single-particle wave functions may be written as 



E*+M* 



JO (>) 



where 



, A E*(E*-M*)R 

4wR 3 j$(uj a )2EfR-2E* + M*' 1 ' ; 



By taking the explicit solutions of the Dirac equation 



E*+M* 



Mrl = N a \ _V^^ |$ femCT , (2.41) 

where the upper (or lower) sign refers to k positive (or negative), it is easily verified that 
only k = 1 (or k — — 1) leads to an angle-independent result on the right-hand side of 
Eq. (2.13). Thus only states with j = 1/2 can satisfy the nonlinear boundary condition 
as given. 



2.3 Radius and Binding Energy of Nuclei 

We have seen that the only change in the calculation of the energy in the MIT bag model 
for nuclei is the addition of a surface term, BS. It is assumed that B is a universal 
constant, chosen to fit one piece of data. Once B is chosen, because of the nonlinear 
boundary condition the radius of the bag is uniquely determined for each nuclei. 

The meaning of this addition to energy-momentum tensor can be clarified by consid- 
ering the total energy of the bag state, 

P° = E(A) = J d 3 x Tg> ag = j d 3 x (T 00 9 v + B6 g ), (2.42) 

which we shall label E(A) as a precursor to our discussion of binding energy later. Based 
on (2.35) and (2.42) the nuclear energy E(A) is derived immediately 

E(A) = AE a + AnR 2 B 



• ,2 



A\l -j^ + M* 2 + AV + 4tcR 2 B. (2.43) 

The first term is the kinetic energy, while the second is a surface term. Essentially it 
implies that it cost an energy BS to make this tension at the bag surface within which 



3 We use tfdr r 2 jl(cor/R) = f + & ±1 (u>) - ^fciH] 
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the baryons move. It should be intuitively clear that energy-momentum conservation is 
related to pressure balance at the bag surface, so that a small change in radius should not 
significantly increase E(A). Nevertheless, the nonlinear boundary condition implies that 

i = 0. (,44, 

We wish to stress that it is an assumption of the model that B should be constant for 
all nuclei. As all nuclear bags have radii in the region (1.0 ~ 1.2) A 1 / 3 fm, this assumption 
will be severely tested. 

Generalizing Eq. (2.43) to include exited states, the nonlinear boundary condition 
implies 



dE(A) 
OR 



Auji 



and hence 



A 



R 2 ^u 2 a + M* 2 R 2 
8nB 



+ 8nRB = 0, 



-R 3 y/u* + M* 2 R 2 



(2.45) 



(2.46) 



The real and positive root R of Eq. (2.46) is found out after an algebraic manipulation, 



R = roA 1 ^, 



where 



1/3 



a = 



P = 



a = 



1/2. 



n — 0, k — —1, 



(2.47) 



(2.48) 



(0/2) 



1/4 



[1 - (/?/2) 3 / 2 ] 1/2 + (/3/2) 3 /^ 



256a \ 



1/2 



27 J 

( AM* 3 ^ 2 



+ 1 



1/3 



/256a\ 



1/2 



1/3 



(2.48) shows that r actually depends weakly on A. 

The above obtained formula (2.47) is well known in nuclear physics. It is one of the 
main successes of our model. 

Using Eq. (2.47) we can then simplify the expression for E(A): 

' \/u 2 a + r 2 M* 2 A 2 / 3 



E(A) = AV + 



r 



+ 4nBr z 



A 2 '\ 



(2.49) 



Hence the binding energy per nucleon is obtained 
e(A) = -(M-0o) + 



A' 1 / 3 . 



(2.50) 



18 



Clearly, the remarkable result obtained for the binding energy per nucleon of the bag 
model of nuclei was indeed a coincidence. 

As was known, the semi-empiric formula of Weizssacker [8, 11] for binding energy per 
nucleon reads 

/ Z — N\ 2 Z 2 
e(A) = -a, + a 2 A^ 3 + a 3 {—^- J + «4^, (2.51) 

in which a±, a 2: 03, and CI4 take the following values, in the energy unit equal to 0.9311 
MeV, 

ai = 16.9177, a 2 = 19.120, a 3 = 101.777, a 4 = 0.7627. 

Now let us indicate that (2.51) is possibly derived from our model if the parameters 
0o ? Vo and B are fitted adequately. Next confronting (2.50) with (2.51) we conclude that 
the above mentioned parameters must fulfil equalities, 



r 



M-V = a u (2.52) 
= a 2 . (2.53) 



It is worth to notice that parameter Vq is explicitly defined by (2.52). The equation 
(2.53) constrains two unknown parameters of the theory, 0o and Vq. 

The final term in Eq. (2.50) represents the contribution to the surface energy from 
the positive-frequency states, where the mass has been shifted by the constant, condensed 
scalar field O - 0o ( or the effective mass M* = M+0 O depends explicitly on the scalar field) 
is a dynamical quantity that must be calculated self-consistently using the thermodynamic 
argument that an isolated system at fixed A and V (and zero temperature) will minimize 
its energy: 

*>>=„, (2.54) 

d(p 

However, O is related to u a by the relation (2.34). As a consequence, (2.54) is replaced 
by conditions: 

[E(A) + \ip] = 0, (2.55) 



duj a 

[E{A) + Ay?] = 0, (2.56) 



in which 



<p = l- M*R - sjul + M* 2 R 2 



tan uj a 

Substituting A from (2.55) and (2.56) one gets finally 

M * R = [M*R+{ul + M* 2 R 2 )]u a t^u a (2 ^ 7) 

uj a tan 2 oo a + (tan oo a - oo a - oo a tan 2 oo a ) ^ool + M* 2 R 2 
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with 0o and u a are the roots of the system (2.34) and (2.57). 

This model problem is exactly solvable. It retains the essential features of the bag 
model for nuclei: A-dependence of nuclear radius and the formula for nuclear binding 
energy. Furthermore, it yields a simple solution to the field equation. This solution 
and model problem thus provide a meaningful starting point for describing the nuclear 
many-body system as well as a consistent basis for considering nuclei with N ^ Z using 
relativistic nuclear physics, the bag model and standard many-body techniques. 

We proceed to investigate nuclei with N ^ Z. 
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Chapter 3 

Bag Model of Nuclei (N ^ Z) 



To realistically discuss nuclei N ^ Z, it is necessary to extend bag model of nuclei N = Z 
to include neutral field, which couple to the isovector current, and the coulomb interaction. 
If Z j£ N, the neutral charged, isovector field corresponding to can develop a classical, 
constant ground-state expectation value in nuclear matter according to 

W> = w'v 

The only change in the results of Chapter 2 is that there are now separate solutions 
for protons and neutrons, with the appropriate frequency modifications: 

Vo + + eA for proton, 



v + U 3 



for neutron. 



The field equation for uniform nuclear matter must also be extend to include contri- 
bution from the classical field 63 and the coulomb potential A : 



il^ - 7V0 - ir 37 % - \e{l + t 3 ) 7 % - (M + O ) 



ip = 0. 



The Lagrangian density, which is obtained from Eq. (2.8) by replacing the classical 
fields. Thus in bag model of nuclei with N 7^ Z 



C = jp 



^^^ fM - 7 V - ^r 37 &3 + ^h°A - (M + O ) 



+ B9 S , 



which is a generalization of Eq. (2.8) to allow for classical, constant fields (f>o, Vo, 63, and 
Aq. Hence, the energy-momentum tensor is still (2.9). 



3.1 Basic Assumptions 

As was know, in Celenza and Shakin theory [2] the motion of nucleon in nuclear matter 
is described by the Dirac equation, 

(^-M-£)V«) = 0, (3.1) 
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where E is the nucleoli self-energy having the form 

£ = o + 7 Vo, (3.2) 

with 0o an d V constants. 

Taking into account the isotopic degree of freedom and the Coulomb interaction of 
protons, in our model, we assume E has the generalized from 



E = 0o + 7°Vb + 7^7% + ^e(l + T 3 ) 7 A), 



(3.3) 



in which O , Vq, &3, and A are constants. Inserting (3.3) into (3.1) we obtain the equation, 

i> = 0, (3.4) 
where the mass matrix M is, of course, 

M = 



^ - 7 °Vb - ir 37 % - ^e(l + t 3 ) 7 A) - (M + 0o) 



M p 



M n 



Now our basic assumption is formulated: the nucleus A is considered to be a M.I.T. 
bag, inside which the motion of nucleon is described by the Dirac equation (3.4). The 
parameters 0o, Vq, 63, and Aq will be fitted to experimental data. 

Let us next consider the energy-momentum conservation for nuclear bag. Let T^ u be 
the energy-momentum tensor of nucleon, described by Eq. (3.4), inside the bag. Then 
the total energy-momentum tensor T-g^ of the nucleus is clearly given by 



T Bag = T " V0V + B9s > 



(3.5) 



where dy and 6s are the well-known step functions for volume and surface of the bag, 
respectively, and B is the surface tension. 

From (3.5) it follows that 



M Bag 



1 



n.d(W)\ s + 2B 



s- 



The energy-momentum conservation, 



leads to 



a = 0, 

^ Bag ' 



B = --n.d(W)\ s , 



(3.6) 



which resembles the nonlinear boundary condition in the bag model for baryon [3]. 
Equation (3.4) and the relation (3.6) constitute the basic ingredients of out model. 
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3.2 A-dependence of Nuclear Radius 



The solution of (3.4) can be found in the form 



^n(x p 



(3.7) 



(3.8) 



Substituting (3.7) and (3.8) into (3.4) we arrive at the equations for protons and 
neutron, separately, 



-ia.V + V Q + -b 3 + eA + (3(M P + O ) 



-ia.V + V - -b 3 + (3{M n + O ) 



E n i> n (x). 



(3.9) 

(3.10) 



For convenience, let us define 



M; = M p + O , 
M* = M„ + O , 

= E p -V -h 3 -eA , 
E* = E n -V + h 3 . 

With this in mind we rewrite (3.9) and (3.10) as follows 

l-ia.V + (3M;]il) p {x) = E* p ^ p (x), 
[-ia.V + (3M*]i; n (x) = E*^ n (x). 



It is known that the solutions of (3.11) and (3.12) read, respectively, 

i/j p (x) = N p 



E*+M* . /Tjr N _ 

3k-i{W p r)$ km 
>\l E * p ~Jf* p jk{W p r)®- km 



N„ 



E*+M* 



3k-i{W n r)$ km 



(3.11) 
(3.12) 



(3.13) 



(3.14) 



where k and k — 1 are the indices of the eigenfunctions corresponding to the eigenvalues 
of operator K, 



K = 0(3.1 + 1), 



w: = e, 



*2 



M* 2 



W. 



K 2 ~ M n, 



(3.15) 
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and the normalization constants N p , N n are defined by 

d z x^ p 9 v = Z, (3.16) 
J d 3 x^J v = N, (3.17) 

Z and N are the numbers of protons and neutrons contained in the nucleus A. 

Now let us assume that the bag has a spherical shape with radius R. Then the 
boundary condition of the M.I.T. bag model is used, which provides the eigenfrequencies 
of proton and neutron, Q p and Q n , correspondingly, 

Q p = W P R, 

Q n = W n R. (3.18) 
As was known, Q p and Q n satisfies the equations 

tanO, = - g =, (3.19) 

tan = =. (3.20) 

1-M*R- y/W n + M* 2 R 2 

It is worth to remember that (3.19) and (3.20) are derived for k — 1, to which corre- 
spond the only states satisfying (3.6). Taking into consideration (3.15) and (3.18) we get 
the energy spectra for proton and neutron, respectively, 



E p = ±\l^ + M; 2 + V + -b 3 + eA , (3.21) 

[Tp i 

En = ±^ + M? + V --b 3 . (3.22) 
For convenience, the sign (— ) drops out in what follows. 

Based on (3.21), (3.22), and (3.5) the nuclear energy E(A) is derived immediately 

Rfi [en Z — N 

E{A) = Z\Jj^ + M; 2 + iV^_| + M* 2 + AV + — ^— h + eZA + AttBR 2 . (3.23) 

Next let us introduce the mean mass (M), the effective mass M* and the mean fre- 
quency (Q) of nucleons contained in nucleus A, 

{M) = ZM p + NM n 



A 

M* = (M) + O , 



A \l^ + M * 2 = Z \l%+ M ; 2 + N \l§ + M n 2 - ( 3 - 24 ) 
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It is easily prove that (Q), defined by (3.24), really exists. 
Substituting of (3.24) into (3.23) leads to 



E(A) = A\I^ + M* 2 + AV + ^-y^&3 + eZA + 4irBR 2 . 

The nonlinear boundary condition (3.6) requires 

dE{A) 



(3.25) 



dR 



0, 



which yields 



„ 871-5 o / fM*R 
A= ^^R 3 \ 1 + 



(3.26) 



The real and positive root R of (3.26) is found out after an algebraic manipulation, 



R = roA 1 /*, 



where 



\4tcB J 

(e/2)V4 



r 

§ l/2 



(3.27) 



(3.28) 



[l_ (e/ 2)3/2]V2 + (e/2) 3/4 



256a 



1/2 

+ 1 ] +1 



1/3 



/256a 



1/2 
11 "I 



1/3 



a = 



/ AM* 3 N 2 
V87r£(ft) 



(3.28) shows that ro actually depends weakly on A. 

The above obtained formula (3.27) is well known in nuclear physics. It is one of the 
main successes of our model. 

Finally, the normalization constant N p and N n given by (3.13) and (3.14), are calcu- 
lated for k — 1, 



I 1/2 



N„ 



( z \ 1/2 [e;{e;-m;)r] 
Wr 3 ) jo (n p ) [2EfR - 2E; + m;] 1/2 ' 

( n \ 1/2 [e:(e:-m:)r] i/2 

WR 3 J j (n n ) \2EfR - 2E* + M*] 1/2 



(3.29) 
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3.3 Weizssacker Formula 



As was known, the semi-empiric formula of Weizssacker [7, 10] for binding energy per 
nucleon reads 

(Z-N) 2 Z 2 

y ' 1 " (3.30) 



/ = -ai + a 2 A- 1/3 + 



a 3 - 



a 4 - 



4A 2 ' ~* A 4/3> 

in which a±, a 2 , a 3 , and a 4 take the following values, in the energy unit equal to 0.9311 
MeV, 

ai = 16.9177, a 2 = 19.120, a 3 = 101.777, 

and 

3e 2 

a 4 = — = 0.7627. 

5r 

The charge distribution radius r for almost nuclei is fitted to be 

r c = 1.2162 10~ 13 cm, 

(3.30) agrees well with experimental data for most nuclei. 

Now let us indicate that (3.30) is possibly derived from our model if the parameters 
0o, V , 63, A , and B are fitted adequately. For this end, let us substitute (3.27) into 
(3.25), 



E(A) = AV + 



r 



A2/3 + b 3 (Z - N) + ^ (3 31) 



Therefrom, the binding energy per nucleon is obtained 
f = V -{M) + 



^{Q} 2 + r%M* 2 A 2 / 3 



,_ 1/3 6s(Z-iV) , eZ 

A 1/3 + V ^ - + A n — . 



■In—- (3.32) 



Next confronting (3.32) with (3.30) we conclude that the above mentioned parameters 
must fulfil equalities 



Vo-(M) = -ai, 



^{Q) 2 + r 2 M* 2 A 2 /* 



+ AitBrl = a 2 , 



b 3 
A 



Z-N 



a 3 - 



2A 
Z 



a 4 



3r eZ 



(3.33) 
(3.34) 

(3.35) 
(3.36) 



It is clear that (3.33) and (3.36) express directly the physical meaning of b 3 and A : 

- The mean field value 63 is proportional to the relative ratio of the numbers of protons 
and neutrons, contained in nucleus A. 
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- The mean value of Coulomb potential created by Z protons equals to that created 
by a sphere of charge Ze, embedded in a nuclear medium, the dielectric coefficient 
of which is 3r /5r c . 



It is worth to notice that three parameters Vq, 63, and Aq are explicitly defined by 
(3.33), (3.35), and (3.36). The equation (3.34) contains two unknown parameters of the 
theory, O an d B. As was shown in the Walecka theory [1], O is a dynamical quality and 
therefore it is defined self-consistently. Namely, we use the thermodynamic argument that 
an isolated system with fixed baryon number A and volume V will minimize its energy, 

*>>=„, (3.37) 

o<Po 

However, O is related to Q p and Q n by the relations (3.19) and (3.20). As a conse- 
quence, (3.37) is replaced by the conditions: 

d 

[E(A) + atiipi + a 2 (f 2 ] = 0, (3.38) 



on 



[E(A) + a 1 <p 1 + a 2 <p 2 \ = 0, (3.39) 



n 



d 

— [E(A) + a 1 ip 1 + a 2 ip 2 ] = 0. (3.40) 



in which 



Vl = i - m;r - ^n 2 p + m; 2 r? - 



i Pl = l-M* n R- + MfR 2 



tan Q p 
tan Q„ 



Eliminating a\ and a 2 from (3.38) and (3.39) and substituting them into (3.40) one 
gets finally 



2ZQ p sm 2 tt p ZM*Ry/Qj + M* 2 R 2 2NQ n sin 2 Q n NM*Ry/Ql + M* 2 R 2 

sin 2 tt p + 2Vt p ~ M;R + y/Ul + M; 2 R 2 + sin 2 fi n + 2Vl n ~ M*R + ^/U 2 n + M* 2 R 2 ~ ' 

(3.41) 

O and Q p , Q n are the roots of the system (3.19), (3.20) and (3.41). 
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Chapter 4 

Conclusion and Discussion 



In the previous sections the basic assumptions and the general results of our bag model 
for nuclei are presented in detail. The model is built on a simple hypothesis: the nucleus 
is considered to be a MIT bag, in which the motion of nucleons is described by the Dirac 
equation and the mean field values <po, Vq, 63, and A$ are supposed to be constants. In 
addition to these mean fields, there exists the surface tension B of the bag that guarantees 
the energy-momentum conservation. 

Two major successes are: the formula (3.27) for the nuclear radius R and the Weizs- 
sacker formula with the parameters verified in (3.33-36). All the parameters appearing in 
the theory <p , Vo, 63, A , and B are, in principle, determined by the equations (3.33-36), 
(3.19, 20), and (3.41). Thus, our theory is a mathematically closed system. The devel- 
opment of the formalism suggested here will be carried out for various concrete nuclei in 
next papers. 
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